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SUMMARY 


During  the  past  several  years,  investigation  has  been  carried 
out  at  the  Daniel  Guggenheim  School  of  Aeronautics,  New  York  University, 
to  study  the  buckling  behavior  of  sandwich  structures.  The  main  effort 
has  been  directed  towards  the  determination  of  the  buckling  loads  of 
circular  sandwich  cylinders  under  different  loading  conditions.  It  was 
found  that  for  cylinders  with  weal:  cores,  linear  theory  predicts  buckling 
loads  which  agree  with  experimental  results.  This  is  contrary  to  the 
case  of  homogeneous  cylinders,  where  linear  theory  gives  buckling  lodds 
much  higher  than  those  observed  from  experiments.  The  reason  for  this 
has  been  explained  in  Reference  11,  In  this  report,  the  theory  of  buck- 
ling of  sandwich  cylinders  und^-v  axial  compression,  torsion,  and  bending, 
and  combined  loads  is  developed  in  a unified  manner;  and  the  principal 
experimental  results  are  presented.  The  interrelationship  obtainod  be- 
tween the  critical  loads  is  plotted  in  the  form  of  non-dimensional  in- 
teraction curves. 
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Symbols  and  Units 
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radius  of  cylinder  to  middle  surface,  in. 

shear  rigidity  of  the  corej  C = (h*t)G  , lb.  per  in, 
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portion  o^  boundary  over  which  displacements  are  prescribed, 
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bending  rigidity  of  a sandwich  element j D ■ 3-t(h*t)  /2(l-<v 
lb. -in. 

Young’s  modulus  of  elasticity,  lb,  per  sq,  tn, 
secant  modulus  of  elasticity,  lb,  per  sq.  in, 

2tE^/a  , lb.  per  cu,  in. 

shear  mod ulus j G = E/2(l+v),  lb.  per  sq.  in. 
thickness  of  core  layer,  in. 

moment  of  inertia  of  sandwich  cylinder  about  its  diameter. 


complementary  energy,  in, -lb, 
buckling  coefficient 

torsional  instability  coefficient 

length  of  cylinder,  in. 
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[l«n */X  (T/7ra2C)  (nx/A  )]  IC/Mat 
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longitudinal  wave  length,  in. 
circumferential  wave  length,  in. 
free  length  of  cylinder,  in. 


m 

M 


wv 

JW*W 


Vy 


number  of  half  waves  longitudinally 
applied  bending  moment,  in, -lb. 

resultant  bending  and  twisting  moments  in  the  composite 
structure,  in. -lb,  per  in. 

x-  and  y-componente  of  the  resultant  moment  on  the  boundary 
in, -lb.  per  in. 

resultant  normal  forces  in  the  composite  structure,  lb. 
per  in. 
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resultant  shearing  force  in  composite  structure,  lb*  per  in* 
resultant  mediar<~ surface  normal  forces,  lb*  per  in. 
resultant  median -surf ace  shearing  force,  lb,  per  in. 
applied  norr>iai  force  in  axial  direction,  lb.  per  in. 
number  of  full  waves  circumferentially 

external  lateral  pressure  on  cylinder  surface  lb,  per  sq.  in* 
l/2(pu-p^)  lb.  per  sq.  in. 

left-hand  side  of  Modified  Donnell* s Equation  with  assumed 
function  v substituted 

shear  stress  resultants  in  composite  structure  normal  to 
middle  surface,  lb.  per  in. 
resultant  Q on  boundary,  lb.  per  in. 

radius  of  curvature  of  median  surface  of  sandwich  structure, 
in. 

bending  stress  ratio,  2Mat/lC 
compression  stress  ratio,  R^51  N-./C 

V -i- 
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torsion  stress  ratio,  R,p*  T/2rra  C 

thickness  of  face  layer,  in, 

applied  torsional  moment  (Fig.  2),  lb. -in. 

displacements  ir.  r:-;  y-,  and  r.-directions  respectively  of 

a point  in  middle  surface  of  cylinder,  in. 

work  done  by  the  force  resultants  due  to  large  deflection, 

in, -lb. 

strain  energy,  inirlb. 

potential  of  external  forces  and  moments,  in. -lb* 
rectangular  coordinates  (Fig.  2) 

x-  and  y-componants  of  the  force  resultant  acting  on  the 
boundary,  lb.  per  in.  1 

curvature  parameter,  2 «*  (l-v^)  1^/2 at 

(h+t)2/a2 

components  of  change  of  slope  of  normal  to  middle  surface 
of  sandwich  shell. 
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Subscripts 


c 

S 

u 

m 

cr 


face  layer 
core  layer 
lower  face  layer 
upper  face  layer 
median  surface 
critical  (buckling) 
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Introdu-’ tion 

A ?ai*d-«d  s-tv-po  structure  consists  of  two  external  layers 
of  thin  high-strength  material  and  a thick  internal  layer  of  light- 
weight material.  The  former  are  usually  called  the  face  layers 
and  the  latter  the  core*  The  basic  advantage  of  a sandwich  struct or* 
is  its  high  bending  rigidity  and  consequently  buckling  strength  and 
its  relatively  light  ./eight.,,  The  use  of  sandw it. h- type  structural 
elements  fci  aircraft  construction  has  received  considerable  amount 
of  attention  with  the  development  of  high-speed  aircraft*  This 
is  because  at  high  speed  the  surface  of  the  plane  must  be  maintained 
smooth  without  buckling  during  flight*  This  is  a requirement  which 
the  ordinary  thicr-sheat-stringer  type  of  construction  cannot,  meet 
without  undue  sacrifice  of  the  weight  economy* 

The  buckling  and  bending  of  sandwich-type  beams  and  flat 
plates  have  beer,  studied  by  many  investigators.  During  the  past 
several  years*  investigations  have  been  carried  cut  at  the  Daniel 
Guggenheim  School  ct  Aeronautics  of  Hew  York  University  to  study 
the  buckling  behavior  of  sandwich  cylinders  under  various  leading 
conditions.  This  repo:’,  summarises  the  principal  theoretical  and 
experimental  findings  of  these  investigations. 
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Tbscretical  Development 


Derivation  of  the  Governing  Equations 

In  order  to  derive  a complete  system  of  equations  for 
the  sandwich  shell  composed  of  face  layers  and  core  layer f it  is 
necessary  to  consider  separately  the  statics  of  the  race  layers 
and  that  of  the  core  layer  of  the  shell,.  Combination  of  the  re- 
suits  obtained  for  the  components  will  lead  to  these  differential 
equations  of  equilibrium  for  the  composite  shello  The  face  layers 
are  treated  like  thin  shells  of  thickness  t having  negligible 
bending  stiffness  about  their  own  middle  surface „ The  loads  applied 
to  these  face  shells  are  of  two  kinds s external  loads  and  loads 
caused  by  the  stresses  in  the  core  layer  (Fig„  l) o The  core  layer, 
of  thickness  h,  is  assumed  to  behave  like  a three-dimensional 
elastic  continuum  in  which  those  stresses  which  are  parallel  to 
the  faces  are  negligible  compared  with  the  transverse  shear  and 
normal  stresses,, 

Assuming  large  deflection,  the  equilibrium  differential 
equations  for  the  upper  face  layer  are  the  followings 
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Similarly*  for  the  lower  face  layer*  the  equilibrium  equa- 
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where  the  subscripts  u ana  i denote  the  upper  and  lower  surfaces  respectively,, 


Under  the  assumption  of  negligible  face-parallel  core  stresses, 
the  equilibrium  equations  for  the  core  layer  are^ 
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Let  the  values  of  the  three  stress  components  at  the  middle  surface  (z=0) 
be  designated  by  the  subscript  su  Integration  of  Eqs.  (7)  to  (9)  results  ill 
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let  and  be  the  resultant  transverse  shear  forces, 
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From  Eqs.  (10).,  (11),  13)  and  (iZ.).  the  following  relation  may  be  obtained*: 
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For  thin  sandwich  shells,  ^ «1.  In  such  cases,  a combination 


of  Eqs.  xl2)„  (13)  r and  (14)  gives 
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In  view  of  the  fact  that  all  face-parallel  core  strej-  & are 
neglected,  the  face -parallel  stress  resultants  and  couples  of  tfa*  composite 
shell  are  due  to  the  stresses  in  the  face^Layers  only  and  may  be  obtainsd 
as  follows: 
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Define  the  deformation  of  the  composite  shell  to  be 
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where  w represents  the  effective  transverse  deflection  of  the  middle  surface 
and  e represents  the  effective  transverse  normal  strain  for  the  composite 
shell.  From  these  definitions „ vu  and  may  be  written  in  terms  of  w and 
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With  equations  (15)  to  (20),  the  following  equations  may  be  obtained 
by  carrying  out  addition!  ps  well  as  subtraction  of  equations  (1)  and  (4), 
and  (^)  and  (5}« 
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When  the  radius  of  curvature  R be cons 2 infinite,  Equations  (21) 
to  (26)  reduce  tc  the  same  equations  obtained  by  Raissner  9 . When  the 
effective  transverse  normal  strain  e for  the  composite  shell  is  neglected 5 
Equations  (21)  to  (25)  reduce  to  the  usual  equations  of  force  and  moment 
equilibrium  for  homogeneous  shells.  Eq.  (26)  has  no  counterpart  in  the 
theory  of  homogeneous  shells,  in  the  sense  that  the  corresponding  equation 
for  the  homogeneous  shell  contains  information  that  ic  not  of  practical 
interest  and  is  therefore  never  formulated.  This  equation  gives  the  local 
change  of  thickness  of  the  shell  caused  directly  by  the  external  loads  by 
way  of  the  non-linear  terms  having  stress  resultants  and  couples  as  factors. 


In  order  to  solve  these  differential  equations  of  equilibriums 
it  is  necessary  to  derive  an  appropriate  system  of  stress-displacement  re- 
lations. This  may  be  done  by  the  use  of  the  method  of  complementary  energy, 
which  states  that  the  true  state  of  stress  is  distinguished  from  all  static- 
ally correct  states  cf  stress  by  the  condition  that  the  complementary  energy 
be  a stationary  value,.  In  the  linear  elasticity  theory,  for  a material 
obeying  Hooke’s  law  and  for  given  surface  stresses  or  displacements  the 
complementary  energy  is  the  difference  of  the  strain  energy  V and  of  the 
virtual  work  W which  the  surface  stresses  do  over  that  portion  of  the  sur- 
face where  the  displacements  are  prescribed.  In  the  non-linear  theory  of 
elasticity  it  can  be  shown^  that  the  expression  for  the  complementary  energy 
may  be  derived  from  the  expression  for  the  potential  energy  by  a Legendre 
type  of  transformation, and  in  the  case  of  a thin  shell  with  finite  deflec- 
tion, the  complementary  energy  J Is 
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J c-  V + D - W (1?) 

where  V is  the  strain  energy,  W,  the  work  function  and  IT,  tho-ifork-dows 
the  stress  resultants  in  the  middle  plane  due  to  large  deflection,  r.-.ji3ly 
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Senbte  the  properties  of  face  layers  and  the  core  by  the  sub- 
scripts f and  c respectively.  The  strain  energy  of  the  face  layers  and 
the  core  may  he  written  as^ 
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The  work  done  by  the  stress  resultant  due  to  large  deflection 
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With  the  aid  of  Equations  (17),  (19)  and  (20),  U becomes 
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The  work  function  W*  is 
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where  the  last  two  integrals  are  the  work  done  by  the  shear  stress  resultant 
in  the  core*  The  first  one  of  these  two  is  the  work  done  due  to  the  dis- 
placement of  the  middle  surface  and  the  second  due  to  compression  of  the. 
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Define  the  diaplaoewuta  of  the  middle  surface  and  the  resultant 
load  on  the  boundary  of  the  composite  shell  to  be  as  follows x 
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V * V * <*n  ? ♦ \ ? 4 » 

Ox  cy 


P*  4 V ^ * (M«  £ 4 V £>' 


-<«nx  f * V ?*>  • 4 ^ <*nf4*nf>4 
ox  oy  4 ox  Oy 


♦ £ii  A ♦ ^)qJ  d . 

6y  j 


12EC  oy 


Introducing  the  Lagraagian  multipliers  ^ to  and  carrying 
out  tho  first  variation  with  respaat  to  N^,  N„,  M^g  Q^-Vnd 


0^  end  az  independently,  the  relation 


I 
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L ■ </  V-  + V + d U + </« 
i c 


9N 

x 

L 9x 


9N  N j,  M , 1 
+ _i2  + JL  „ J£  £e  ! 

9y  R 9x  R 6x  J 


rSN  9N  N 

_J2  + *JL  _ _ 

9x  9y 


jsz  fist  + !!&:  2s.  _ if] 

R 9x  R 9x  RJ 


1 [Hi  ,L(,  S!),i  (H 

1r  3x  x9x  9y  ^ 


9v^ 

9x 


♦-(«-)+-  (K  - (MT  24)  - 2-  (m  2») 

9x  ^ oy  dy  y dy  9:c  ' 9x  9y  ^ 9x 


9Q  ''<1 

- •*-  (M  % . 2_  („  SS)  , p + (_2  , li) 

9x  ^ 9y  9y  y 9y  9x  9y 


+ (q.  ~ + q ^ i & 


<V~) 

9x  y 5y  R 9yJ 


„ T9M  9M  M _ ,.,.*2  N 

+ 1 _x  + _JZ  + _Z  9^  _ IkllL  _Z  9e 

' i_  Ox  9y  R 9x  4 


R 9x 


X 


0M 


*Z  + fUz  _ 1j2:  9w  + liLtil2.  .JZ  .qe  __ 


j_  9x  9y  R 9x 


R 9x 


X 


M 


+ - (M„  2s  ♦ M.„  2S)  t 2.  2s  + M &) 


L R 9x  9x 


v 9y  9y  ^ 9x 


9y 


i i 


(hvt)2  a 


; M --- 

•V 

av  ■*  ax 


N 2_  (N 

*y  *y 


— + N 22.) 


~ /,  . \ 2 0 

/,.>>'  v / „ dw  . „ 0W\  . y 

\h+S.;  \o..tq;  - (Q. . " • v Q._  — J + ■>-.  — « — —*• 
x ax  * dy  4 R ■ 


dx  dy  r-  0 


Lead'.-  to  the  conditions 


X,  - 


X2  - " 


A 3 “ v 


^4  “ 


^5  = py 


X6  - e 


c«n  the  boundary.  As  equation  (32)  also  holds  for  any  part  of  the  structure 
if  the  boundary  displacements  referring  to  this  part  Ere  identified  with 
the  displacement  occurring  in  the  actual  solution  of  the  problem,  it  fellows 
that  the  Lagrange's  multipliers  throughout  the  structure  are  related  to  the 
generalized  displacements  in  the  interior  of  the  boundary.  It  folxovs  that 
equation  (dd)  also  holds  true  now  in  the  interior. 

Using  equations  (33),  the  Euler  equations  of  L = 0 in  the  Cal- 
culus of  Variations  give  the  following  stress-displacement  relations: 


1 [~  a- . 

-i.  i ' 'J  w \ 

“ j « — ) 


I ' 

3x  .c  L dx 


th+i)2  rdex2  i 
y—) 

4 dx  J 


9Z  _ 2 + 1 ["(§2)  + (h+t)2  ^6e^ 


dy  R 2 I dy 


u 9v 
R 3x 


de 


R dx 
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2{i+<j)N^r  ..  3u  + 3r+9wdv  + (h+t)2  6e  8e 
2tE  0y  dx  dx  3 y 4 8x  3y 

+ v 3w  (h+t)2  ^y  3a 
R 8x  4 R 8x 


M - vji 

x i y .... 

^x 

B 

+ 

8e 

t(h+t)^Ex/2 

dx 

3x 

8x 

My  - VA 

+ &L 

* 

u 

8s_fx 

SB* 

+ s. 

— 

t(h-M)2Ef/2 

0y 

3y 

3y 

R 

dx  R 

3x 

R 

4(l+v,)M 

.!5*H 

h !5s 

dw 

& 

dw 

S« 

V 

i & 

t,.. 

6y 

3x 

ox 

8y 

8y 

8x 

R 

8x 

R 3x 

Sl _ Ikikl  k.  (fi  + £i)  + 9w_3«w_8£  e 

(h+t)Gc  12E(>  8x  5x  oy  X dx  dx  ~ dx 


Sz 

(»t)04 


h+t 

3J2E 

c 


+ = <?+§£ 
3y  8x  8y  ^ 5y 


8e 

— v 

9y 


v + v dw  + (h+t)4  e 

R R 8y  4 R 


e « '•  <j 


/E 

znr  u 


(34) 


In  the  case  of  flat  sandwich  plates*  R—*"00,  the  stress'-strain 

2 3 

relations  (34)  reduce  to  exactly  those  obtained  by  Reissner  '' . In  the 
case  of  homogeneous  shells  (G  « “ E — ► °o,  e-*0),  equations  (34)  then 

C «t  C 
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rodace  to  the  stress-strain  relations  for  homogeneous  shells  when  finite 
deflection  is  assumed,, 

If  buckling  of  the  sandwich  structure  as  a unit  is  to  be  studied.,, 
the  effect  of  the  transverse  normal  stres»  d^fdmsabili'by  6aji  usually 
neglected  and  therefore  one  may  assume  Eq  * oo  and  e = 0o  In  the  buckling 
theory,  the  displacements  u,  v,  and  v are  the  displacements  occurring  during 

buckling  and  the  quantities  N 9 N , etc.  represent  only  the  changes  in  the 

' ^ 7 

internal  forces  during  buckling.  If  we  denote  the  median  surface  stresses 

- 5 

by  N 9 N « and  N and  neglect  the  terms  which  Donnell  ' justified  as  being 
x y jty 

small,  the  equilibrium  equations  beoomB 


3N  9N 

x , xy  ^ 

“ + R 0 

dx  3y 


9N  3N 

— TZ  o 

dx  3 y 


!i  + 5:_5 ci!r'+;  & 

dx.  dy  R x dx* 


— 3^w  *"  Q^w 


(35) 

(36) 


(37) 


3M  3M 

~ + ~ 22  - Q_  » 0 (3B) 

dx  dy 

0M  3M„ 

+ - Q,.  c 0 (39) 

dx  0y  ^ 


and  the  stress-displacement  relations  become 


N - v JJ  - 

jjt x_y  an 


2tE. 


dx 


2t,E, 


A 


i dv  w 

R 30  R 


(40) 
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—<-V  — 


2(1+VV)N 

ZLjsl 
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1 0u  d? 

R ac  3x 


M - v M 
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t(h+t)2Ef//2  9x 


a 6. 


^ 4C  c o nt* 


M — v 
Jl _LJL 


as 


4(l+vf)M. 
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t(h+t)2E 


!5s  .1 


07 


ao 


«* 


(h+t)G 


= 8 + 
^x 


dw 


ox 


2t 


<h+t)Oc 


1 

« 6 ■»  -•  — 
y r as 


Equation  (AO)  may  be  rewritten  in  the  following  fora: 


N » » 

x i 2 

l"Vj. 


2»  * v,.  (A  & 

dx  R 94 


W)1 

R J 


N 


2t.E  - * 

r ^ 

f 


0U 

3x 


+ i & _ w 1 
R 0fi  R J 


2tE, 


N 


/ 1 9u  dv 


^ 2(l+vf)  R 60  8x 


■) 


(41) 
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(h+t.)  tEr  0|3 


x 2- 

2(l-vf  j 


x 1 v 

l — — + v , *■  ; 

'ax  r as 
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v 


(h+t)*tE*  S?.r  , 

± t X L 1 \ 


»PL 


7 2(1  -o 


9x  R 3D 


(h+t)2t£,  . as  9R 

M k ( - J*  + .11.) 


*7 


2(1  + vf)  R lJ0  dx 


Q*  ■ (h+t)  G,;  (px  + 


8y  % 
0x 


Q_-  - (h+t)  G (p  + * “) 

C 5 R dd 


(41  oont’dc) 


Substituting  these  relations  into  equations  (35)  through  (39) 
and  eliminating  a.  v.  B„.  ar.d  p,P  one  obtains  finally 

-V 


D V°v  + (1  - - V2) 


2t£i  A 


L R 


?>x 


- V4(HX^^2K  ^ vN  4;}1  - 

3x  v oxdy  oy  J 


0 


(42) 


where  V2  - QZ/dx*  + 02/dy2,  D « Eft(h+t)2/2(l  - vf2) 
and  C ■ (h+t)  Gr, 

Equation  (42)  is  the  Modified  Donnell *3  Equation  for  curved  sandwich  plates 
ana  sandwich  shells  and  is  the  one  which  was  obtairied  by  Stein  and  Mayers^ 
for  ail  isotropic  sandwich  element. 

Galerkin’s  Method  of  Approximate  Solutions 

We  shall  now  proceed  to  solve  equation  (42)  in  the  case  of  buckling 
of  sandwich  circular  cylinders  under  i.'arious  loadings.  Let  us  consider 
the  cylinder  to  be  under  ths  general  combined  loading  of  axial  compression, 
torsion  and  bending  (Fig«  2).  In  such  a case,  we  have 


i 
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N *e  _ JIj  + (2Mat/l } cos  0, 

(43) 

N = T/^rra^,,  N " 0 

*y  f y 

where  is  the  force  per  unit  length  due  to  axial  compression,  M is  tbs 
bending  moment,  I is  the  moment  of  inertia  of  the  cylinder,  T io  th» 
torsional  moment,  and  a the  radius  of the  circular  section.  With  such  a 
loading,  equation  (42)  thus  becomes 


D V8v  + (1  - - V2) 
C 


l 


57 


V 


( - N-  + 


2Mat  „og  afw  + T_  82v 
I 8x2  TTe?  9x90 


] 

>J) 


(44) 

where  y is  i-sJksn  as  aQ* 

Equation  (44)  may  be  solved  by  means  of  Galerkin's  method  as 
follows*  we  first  assume  the  deflection  w of  the  cylinder  after  buckling 
in  the  form  of  a series  that  satisfies  the  boundary  conditions  but  with 
undetermined  parameters.  For  long  cylinders,  nowever.  the  boundary  condi- 
tions at  the  two  ends  become  unimportant  and  we  may  assume  the  deflection 
series  without  any  regard  for  the  end  conditions.  For  a cylinder  under 
combined  axial  compression  and  bending,  w may  be  assumed  in  the  following 

,___7 

XULiU  f 

v ■ sin  £ A cos  nO  (45) 

1 nO 


and  for  a cylinder  under  torsion  only,  the  deflection  after  buckling  is  of 
the  following  form  , 


u 


oO 

T7 

XJ 

n«l 

oo 


b sir  ^ - n£>) 
n 1 


(46) 


£ B (sin  coo  nG 

n ^ 


cos  sin  rO) 
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Guided  by  these  expressions,  ve  shell  assume  the  deflection  of  the  cylinder 
in  the  case  of  combined  loading  as  follows; 


sin 


m7Tx: 


n- 


X 

n**0 


A cob  + ccs 
n 


rarx 


nl 

£ 

n^l 


E sin  nQ. 
n 


(47) 


With  v expressed  in  a proper  series,  vs  shall  next  substitute  this  series 
into  (44)®  If  the  expression  (47)  happens  to  be  the  exact  solution  of  equa- 
tion (44),  after  substitution  equation  (44)  will  be  identically  equal  to 
zeroo  In  general^  this  will  not  be  so  and  the  resulting  expression  will 
be  a function  of  x and  0 which  we  shall  denote  by  Q.  Galerkin:s  equations 
for  the  determination  of  the  coefficients  A.  and  B are 


r 


Q sin  cos  r.O  ad©  dx  * 0 


Q co3  sin  nO  adO  dx  = 0 

* 


(48) 


(49) 


which,  when  written  out,  becomes 


tt  9^w  DF  r-,2  d^v\  . mrrx  n . 

F —r  - — V — 7)  cos  nO  adO  dx 

ax4  c ax4  i 


D 6 A2v.  mrrx 

- - v ain 
c dx  s 


cos  n©  £d©  dx 


(50) 
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2 2 

(v4  - “ V6  “-?)  siE  ” COB  n0  cos  0 adO  dx: 


, 4 9v  D 6 Q^w  s mFv  _ jo  j n 

( ^ “ y — — / sin  ~ — ■ cos  nC  adO  dx  » 0P 

3x00  C 3x30  i. 


/ ; 

S r\  < 


in  which  F *=  2tE^/a  , 


2rr  \ 

f / ( D V8w  + ^ “"T  - ^ V*-  ”r)  003  sin  adO  3* 

J j 3x4  C 6x^  V 


0 0 


/ 4 3*w  D __6  8 u,  nrrx 

( v — p - " V — p)  c°s  

3x  C 3x  i, 


?in  nO  adO  dx 


(L  32v  D 6 3^v»  mrx 

<V  “2  “ - V “ 2>  C0S 

3x  C 3x  \ 


sin  nO  cos  C adC  dx 


/ /_4  02w  D 6 82w  n 

/ (V  ' — - “ V — — ) cos 

' 0x30  C 3x" 


sin  nO  ad©  dx  **  0 


(51) 


Substituting  the  assumed  function  v into  these  equations  and 
carrying  out  the  integration,  we  obtain 


i 


.nstead  of  working  with  these  two  equations,  it  is  found  conven- 
ient to  corabine  them  into  the  following  ones*  Thus,  by  subtracting  equation 
(53)  from  (5£),  we  obtain 


I 


(54) 

in  which  ^ = wh'r7A  , Kq  n - E^,  and  ir.  which  2Kq  is  to  bo  substituted 
for  Kr  at  n » 1}  K ■ 0 when  n>  n(9  When  equations  (52)  and  (53)  are 
added,  we  have 


(55) 


i 
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in  which  X*  * A ♦ B Again  at  n = i«  2K’.-.  is  to  be  substituted  for  K1 
nun  u 

and  K!  - 0 when  n > n,  3 
n .1 

The  substitution  of  n from  0 to  in  equations  (54)  and  (55) 
results  in  l^n^  + 1 simultaneous  algebraic  equations  for  2n^  + 1 unknowns,. 
One  solution  to  these  equations  is?of  course,  the  trivial  one,  namely, 

Kn  » K’n  ■ C (n  = 0,  1,  2,  3,  « o o,  n1)D 


Ths  non-trivial  solution  is  carried  out  in  the  following  section 


I 


~2S= 


sion  of  the  Buckling  Load? 


Caso  1.  Buckling  under  axial  compression. 

If  the  sandwich  cylinder  is  under  the  action  of  axial  compress!; 
only,  then  M - T = 0 in  Eqs.  (5!)  and  (55)  and  we  obtain 


D(A“+r7r  . FX4  . DFA4(A  +n  ) 

§ T 4 76 

a a Ca 


2 M2<A2+*2)2 


2,  ->2.  2x2  -1 


n7dA  ( A' V)‘ 
T 


Ca 


“n  0 


(56) 


To  find  the  non- trivial  solution?  of  the  above  equation,  we  roast  set  the 

coefficient  of  equal  to  zero,  which  gives,  after  some  s impl if ic ati on  D 

n2  i 


N1  Ta(l+  ^2)4  +4  T(l-0  ^4  + 2a  \2  (1+  ^2) 

C 2^(l-v  2)  1_(1+  s if  + Ta  (1+  r,2)5 

-*•  r\  **.  ~ 


A2  A 


A" 


(57) 


9 t •j 

where  T r-  (h+t)G,/E^t,  a «•-  (hvt)  /a"  . 


In  view  of  the  presence  cf  the  sandwich  cylinder  parameters  IT 
and  a in  equation  (57),  the  determination  of  a mlrixrum  with  respect 
to  the  parameters  A and  n proves  to  be  intractable  analytically.  Con~ 
sequently,  a numerical  treatment  was  attempted.  In  the  numerical  analysis 
three  cases  were  treated?  (a)  the  acre  considered  extremely  weak  in.  shear 
with  a = 10  and  V «*  2.7  x 10  , (b)  the  core  considered  moderately  week 

in  shear  with  a c 1C  4 and  I ^ 2.7  x 10“  $ and  (c)  the  core  considered 
relatively  strong  in  shear  with  a •-  10"“*  and  If 2.7  x 10“  . 

Tne  values  In  caso  [a)  were  chosen  so  as  '>0  coincide  with  the 
nominal  properties  of  the  cylinders  used  in  the  experimental  study*  The 
other  values  were  chosen  as  convenient,  multiples  of  tr.e  first  set*  The 
results  of  the  numerical  analysis  are  shown  in  Figs.  3,4,  and  5.  The  values 
cf  r.  and  X obtained  from  these  figures  for  a minimum  K,/C  are  given  in 
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Table  1. 


Table  1 


c 

T 

Ya/ir 

N-/C  . 
r rain 

T- 

A 

i<r4> 

2.7xl0=^ 

3.7 

1.00 

o 

10"4 

2.7x10  ~ 

0.37 

0.626 

0 

18 

10~4 

2,7x1a"1 

0.037 

0,0769 

7 

7 

In  Figs.  3 to  5,  curves.  for  only  a few  values  of  n are  shown 
which  are  sufficient  to  ensure  a minimum  value  of  N-j/C  for  the  corre- 
sponding A • It  was  possible  to  rake  approximate  analytical  calculations 
that  served  as  a guide  to  bracket  the  minimum  n for  a particular  X • 

This  reduced  the  amount,  of  calculation  to  rhe  point-  where  it  was  only 
necessary  to  calculate  from  Eq.  (57),  the  n curve  for  a corresponding  \ 
range.  With  this  procedure,  the  values  of  N-,/C^_  are  reasonably  cer- 
tain, although  r.  and  X may  be  ir-  ."Light  error. 


These  numerical  calculations  indicated  that  the  minimum  values 
of  occurred  at  n^  0 andX^Xl®  By  taking  n = 0 in  (57),  we  have 

4-  2aA2+4T(i-vJ) 

(58) 


Ta  \ + 2a  X +4|tl-',j.) 
C 


raX4*2r2(l-vj!A: 


Since  minimum  N-,  occurs  at 
mum  from  the  condition 


large  value 


au/a  A2  » o 


we  can  determine  the  mini- 

( TO'i 

\ + ' i 


Equation  (59)  has  the  following  real  solutions: 


X - °° 


f - 

^ & 1 - OL 

■r'd-v2) 


(60) 

(61) 


Note  that  from  eq,  (6i)  when 
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V 

G a 
c 


i-V 


(62) 


the  only  real  root  is  - oo  . 

The  normal  stress  in  the  faces  at  buckling  is 

o = N,/2t 
cr  1 

Substituting  Eqs.  (61)  into  Eq.  (5?)  and  simplifying,  we  obtain 

m K*r[  (h+t)/al  (63) 

WA  M X i.  J 


where  the  buckling  coefficient  kg  is 


+ _JL 
^2 , 


Eft 


4(l-vf<)a  X"(h+t)  2(l-vf2)Gca 

Sft  ^2  (h+t) 


(64) 


2(l-v.,2)Goa 


Values  for  the  buckling  coefficient  kg  and  half  wave  length  have 
been  computed  from  Eqs,  (64)  and  (61),  respectively,  by  taking  = 0,3, 

These  data  are  shown  in  Figs*  6 end  7*  When  E^t/f  a ^ 0*95,  the  value 
of  ohe  buckling  coefficient  is 

ks  *=  (1/2)  (Gca/Eft)  (65) 

Under  these  conditions,  the  buckling  stress 

0cr  » (h+t)Gc/2t  (66) 

has  a value  that  is  independent  of  the  wave  length. 

An  interesting  feature  of  the  theories  illustrated  in  Fig,  6 
is  the  fact  there  are  apparently  three  ranges  into  which  the  shear  effects 
of  the  core  can  be  divided:  (a)  When  E^t/G^a^.0,95,  the  core  is  extremely 
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weak  in  shear  and  Ncr  ■ (h+t)Gc«  (b)  When  0.1  E^tA^a  ^ 0.95,  the  core 
is  moderatley  weak  in  shear  and  acr  is  given  by  Eq.  (60).  (c)  When  E^t/G  a 
< 0.1,  the  core  is  relatively  strong  in  shear,  and  the  sandwich  cylinder 
behaves  essentially  in  the  same  manner  as  the  homogeneous  cylinder.  While 
shear  effects  arc  negligible  in  this  region,  tire  nonlinear  theory,  which  io 
necessary  to  describe  buckling  of  a homogeneous  cylinder,  will  also  be  re- 
quired  for  the  sandwich  cylinder.  The  linearized  solution  for  the  buckling 
stress  of  the  latter  can  be  obtained  as 

o - (]/\/l-v2)E(h+t)/a 

The  strain  energy  solution  obtained  by  Leggett  and  Hopkins'*'0  also  agrees 
with  the  results  shown  in  Fig.  6 when  their  solution  is  simplified  for  a 
core  that  does  not  carry  any  axis!  load.  A slight  difficulty  in  their 
solution  which  was  not  noted  by  Leggett  and  Hopkins  is  that,  in  the  region 
A-A'  shown  in  Fig.  6,  the  buckling  stress  given,  by  their  equation  goes  to 
zero,  for  in  this  region  X becomes  imaginary  for  this  branch  and,  conse- 
quently, the  X **  °°  branch  must  be  used. 


Case  2.  Buckling  under  torsion 


If  the  cylinder  is  under  the  action  of  torsional  loads  only, 
then  in  Eqs.  (54)  and  (55)  and  we  obtain 
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(67) 


To  find  the  non-trivial  solutions,  we  set  the  coefficient  of  K equal  to 

n 

zero  and  obtain 


+ a T *2CX.  A ^(1-*^-) 

r2A2  t(i4)2*  2*X  * 4 ?(lt  fe)3 

A A AX 


2 

a o 


(68) 


where  the  plus  and  minus  signs  indicate  that  the  direction  of  T is  immaterial 
to  the  buckling, 

In  order  to  find  the  lowest  possible  value  of  T,  we  may  minimize 
the  right-hand  side  of  (68)  with  respect  to  the  variables  X and  (n/  X)  • 

It  may  be  pointed  out  that  X and  (n/  X ) actually  should  take  on  only  dis- 
crete values  on  account  of  the  integral  character  of  the  n and  m 
However,  in  the  case  of  sandwich  cylinders  wit:,  weak  cores,  it  will  be 
shown  subsequently  that  the  minimum  values  of  T occur  at  very  large  values 
of  X . It  is  therefore  possible  to  consider  both  X and  (n/A  ) as  con- 
tinuous variables  in  the  minimizing  process. 

Let  us  first  examine  the  condition  9T/3A"  0.  We  find  from 
this  condition  that  either  X ® oo  or 


2(i-v2f)T  Va  - y^f  T (1+  h) 


Vocd+sJ) 

X* 


o o 22 
(1-v2)  r2(l+  h)  -a 

A 


It  is  evident  that  when 


r 


>Vl-v2f  (l*h) 


the  only  real  root  is  X ■ oo  • 

The  condition  (70)  corresponds  to  the  case  where  the  sandwich 
has  a core  weak  in  shear.  By  substituting  X ■ in  equation  (68)  and 


minimizing  T with  respect  to  (n/ A)»  find 


which  gives 


r T . x 

k 2 ' 

2r,*C  cr 


-2E  - 1 i 


! 
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This  result  fo-»-  a core  very  weak  in  shear  corresponds  to  that  obtained 
for  a sandwich  cylinder  under  axial  compression*  a sandwich  plate  under 
compression,  or  a column.  The  critical  loading  in  all  these  cases  depends 
solely  upon  the  shear  rigidity  of  tbs  cores 

For  sandvich  cylinders  with  cores  not  weak  in  shear,  it  is  neces- 
sary to  proceed  directly  with  the  minimization  of  the  Eq»  (oa)  with  respect 
to  A and  the  ratio  (n/  A )•  This  results  in  a pair  of  non-linear  equations 
in  A and  (n/  A ) and  their  solutions  become  very  tedious.  Inasmuch  as 
it  is  known  that  ths  linear  theory  does  not  predict  adequately  the  buckling 
ioad3  for  a homogeneous  cylinder,  it  is  expected  that  the  linear  theory, 
on  which  the  present  analysis  is  basedj  also  will  not  predict  the  buckling 
loads  for  sandwich  cylinders  with  cores  strong  in  shear,  Hie  reason  why 
linear  theory  checks  well  with  the  experimental  results  in  the  case  of 
sandwich  cylindea*  with  weak  core  has  been  presented  in  Reference  11* 


Case  3,  Buckling  under  bending 

If  the  cylinder  is  under  the  action  of  bending  loads  only,  V9 


let  N^*T=0  in  Eqs,  (54)  and  (55)  and  we  obtain 


D(  X2*°2)4  + F A 4 DF  ).4;  A2+n2) 
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0 


(73) 


for  n f 1 where  0.  For  n 8 1,  2Kq  is  to  be  substituted  for  end 

for  n > K ■ 0,  For  any  value  of  taken  we  then  have  (n^+l)  equa- 

tions with  (n^+l)  unknowns.  To  obtain  the  non-trivial  solution,  wo  shall 

set  the  determinant  of  the  coefficients  cf  K„.  K, . K^.  . . . K , eounl  to 

u’  a-  7 nj 

zero,  from  which  the  buckling  load  can  be  calculated.  This  procedure, 
however,  requires  very  tedious  calculation  for  a general  sandwich  cylinder* 
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To  simplify  the  calculations,  let  us  study  the  effect  of  A 
on  the  buckling  load  of  a sandwich  cylinder  with  weak  core.  This  is  done 
by  calculating  the  buckling  icad  or  such  a cylinder  with  n-^  = 1,2, 3,4, 5. 
The  results  are  clotted  in  Fig.  8.  It  is  evident  from  these  curves  that 
for  all  those  values  of  n^,  minimum  M occurs  at  A- * °o.  The  convergence 
of  3D/2Mat  as  n^  increases  can  easily  be  seen  from  the  figure. 

By  assuming  that  A is  a Isrge  number,  many  terms  in  eq.  (73) 
can  be  neglected  and  we  obtain 


n-1 


+ L 


+ K 


n+1 


(74) 


for  n i 1.  L « Cl/!!at  in  tlie  above  equation.  For  n ■ 0,  we  have 


I*0  + « 0 (75) 

and,  for  n = 1,  we  have 

2KQ  ♦ I*,  + K2  = 0 (76) 

To  calculate  the  buckling  load,  instead  of  solving  the  determinant  equa- 
tion for  a large  number  of  terms,  an  alternative  solution  can  be  obtained 
as  follows.  We  observe  that  Eq.  (74)  is  a finite  difference  equation  in 
terms  of  Kn«  The  trivial  solution  of  Eqs.  (74),  (75)  and  (76)  is 

Kn  = C (n  B 0,1, 2, 3,  . . . ) 

The  nontrivial  solution  of  Eq.  (74)  is 


Kn  - Aan  + Bcf"  (n  = 1,2,3,  . . . ) (77) 

where  A and  B are  arbitrary  constants  and  a is  related  to  L by  the  formula 

+ dL  + 1 c 0 (78) 

or 

I*  " - (a  + a"1)  (79) 

Substituting  Eq.  (77)  into  Eq.  (76)  and  using  relation  (79), we  fiud 
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2K0  - - US.-K2 


- L (Aa  + Ba”1)  - (Aa2  + Bcf2) 


- -A  (a2  + aL)  - BCa***  + a-xL) 


- A + B (80) 

Now,  equation  (75)  is  the  only  one  yet  to  be  satisfied.  Substituting  Kq 
given  by  Eq.  (80)  and  X-  given  by  Eq.  (77)  into  Eq.  (75)  and  using  rela- 
tion (79) » we  obtain 


(A  - B)  (a  - a”1)  = 0 

for  which  we  find  a2  *»  1 and  a ~ - 1.  Hence,  L = - 2 or 


2Mat/l  = 1 C = - (h+t)Gc  (81) 

We  note  that,  when  o = - 1,  the  two  roots  of  Eq.  (78)  are  equal  and  Eq. 

(77)  must  assume  the  form 

Kq  ■ (A  + nB)an 

But  must  be  finite  when  n— f we  have,  therefore,  B c 0,  In  this 
case,  we  find  that  a - - 1 is  again  the  correct  solution.  Hie  plus  or 
minus  3igns  in  Eq.  (81)  indicate  that  a positive  or  negative  moment  on  the 
cylinder  will  make  no  difference  so  far  as  buckling  is  concerned.  Eq.  (81) 
indicates  that  the  sandwich  cylinder  buckles  when  the  maximum  tending  stre3 
reaches  the  buckling  stress  when  the  cylinder  is  under  axial  compression 
only. 


Case  4.  Buckling  under  combined  axial  compression  and  bending 


The  determination  of  the  buckling  load  for  a sandwich  cylinder 
under  combined  compression  and  bending  can  be  carried  out  in  a manner  simi- 
lar- to  that  employed  in  the  case  of  buckling  under  pure  bending.  Let 
T *»  0 in  Eqs.  (54)  and  (55).  Again,  by  assuming  large  X » Eqs.  (54)  and 
(55)  become 


K . + 
n~l 


1 

L K_ 


n 


K - » 0 
n+1 


(82) 
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for  n / 1.  where  L*  = n i - IC/l'at-  For  n * 0,  we  have 

A0  + K1  “ 0 (83) 

and,  for  n = 1,  we  have 

2Kq  + lh1  + K2  * 0 (84) 

Eqs,  (82),  (83);  and  (84)  are  exactly  the  same  as  Eqs.  (74),  (75),  and 
(76).  The  solution  is  therefore  ■=  - 2 or 

Nx  - (2Mat/l)  « C - (h+t)Gc  (85) 

Equation  (85)  again  indicates  that  buckling  of  the  sandwich  cylinder  will 
occur  when  the  maximum  compressive  stress  in  the  cylinder  is  equal  to  the 
buckling  stress  when  the  cylinder  is  under  compression  alcne. 


Case  5.  Buckling  under  combined  axial  compression  and  torsion 


When  M «=  0 and  again  assuming  that  A is  large  while  the  ratio 
n/X  may  be  finite,  Eqs.  (54)  and  (55)  become 
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(86) 


and 


2 Ni  <p 

a*  ^ + 1 


S ) K * *=  0 
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The  non -trivial  solution  occurs  when 

T r.  + (.  nf 

2„  i ’ 1 -\2  “ ) 

Ta  0 A X c 

Minimizing  T or  with  respect  to  the  ratio  (n/X  ),  we  find 


(88) 


n b + T 
X 2ira2C 
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and  buckling  occurs  when 


2tte*T! 
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2rra2C 
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Case  6,  Buckling  under  combined  axial  compression,  bending  and 
torsion. 

Now  let  us  take  up  the  general  case  of  buckling  under  combined 
compression,  bending  and  torsion.  Our  problem  is  then  to  find  the  non- 
trivial solutions  of  Eqs.  (54)  and  (55).  This  is,  however,  a very  difficult 
task.  From  the  previous  cases  of  buckling  of  sandwich  cylinders  with  weak 
core  under  compression  alone,  bending  alone,  and  torsion  alone,  it  was  found 
that  the  minimum  buckling  loads  occur  at  A * oo.  Dividing  equations  (54) 

g 

and  (55)  by  A , and  remembering  that  n^  may  also  be  very  large,  we  shall  drop 
all  terms  containing  X to  a power  greater  than  zero  in  the  denominator 
but  keep  terms  containing  the  ratio  n-j/A  • Equations  (54)  and  (55)  thus 
become 

(14-V^o 

(1  - i-  be  a>  * “£ (K  - * * w ■ 0 <92) 

or, dividing  these  equations  through  by  Mat/'lC,  vs  obtain 
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in  which 
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JLa  in  cases  3 and  4»  one  may  attempt  to  solve  equations  (93)  and  (94)  as 
finite-difference  equations.  This,  however,  was  not  found  possibles  The 
reason  is  that  equations  (93)  arid  (94)  are  now  finite-difference  equations 
with  variable  coefficients  and  the  solutions  of  such  equations  are  difficult 
mathematical  problems.  Instead  we  shall  solve  the  problem  in  the  following 
manner. 

Let  us  first  investigate  the  magnitude  of  n^.  If  n^  is  small 
compared  to  A then  n^/A  — C as  A °°  and  equations  (93)  and  (94) 
reduce  to  the  governing  equations  in  the  case  of  combined  compression  and 
bending.  In  order  that  the  torsional  load  may  have  any  effect  on  the  buck- 
ling phenomena,  must  'also  be  a large  number  so  that  the  ratio  n^/A 
may  remain  finite.  For  any  n close  to  n^,  say,  n = n^-1,  n^-2,  . . ., 
n--q,  (q  n^),  it  is  obvious  that  such  terms  as  l/A  , 2/ A , • . . 

q/  A vanish  as  A — " ► oq.  Thus,  for  the  equations  in  which  n «=  u-,  , n^-1, 
n.,-2,  . . . .,  n^-q,  the  following  relations  hold: 


\ ■ Lnrl  * \-2  * * • * ‘ 
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Therefore,  if  ir,  the  series  for  v-  we  take  the  summation  of  terms 
from  t.n  n^;  the  -system  of  equations  will  be  as  follows? 

UL  ■*  K - *=  0 
nl  nl“1 


K + IX  , + K 0 » 0 
n^-1  n^-2 


K ,.  +&  - 0 

ru-a-rl  n-^— q 
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and 
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where  L'  ■ (1  + ^ 
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To  obtain  the  non-trivial  solution  of  (98)  we  set  the  determinant 
of  the  coefficients  Kn  equal  to  zero,  namely, 


L 1 0 

1 L 1 
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Equation  (100)  can  be  rewritten  as 
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(100) 
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By  tehing  moro  and  mors  heimS  in  the  expression  for  w,  we  have 
successively  larger  and  larger  c >ntinued  fractions  in  (101)  aixi  in  each 
case  the  roots  L are  found. 
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The  question  now  arises  as  to  which  of  the  roots  obtained  in 
solving  equation  (101)  is  the  one  which  should  be  used  in  the  final  anal- 
ysis. This  is  easily  determinable  from  relation  (??),  which,  when  re- 
written, becomes 
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— n.  N 
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Obviously  the  lowest  value  of  the  critical  bending  moment  M occurs  for 
the  largest  value  of  Lj  therefore,  in  solving  equation(lOl) , the  largest 
root  obtained  in  each  case  is  the  one  selected. 

It  is  noted  (Fig.  9)  that  the  largest  root  of  equation  (100) 
approaches  the  value  of  2.00  as  larger  and  larger  continued  fractions  ere 
considered.  Therefore  the  value  L = 2.00  is  taken  as  the  solution.  Equa- 
tion (102)  therefore  becomes,  after  dividing  through  by  Ic/Latf 
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In  order  to  find  the  lowest  possible  value  of  the  term  2Mat/lC, 
it  is  necessary  to  minimize  the  right-hand  side  ">f  (103)  with  respect  to 
the  variable  ri^/\  . Thus,  we  obtain 
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Hence 
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Consequently,  equation  (103)  becomes 
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Exactly  the  same  result  Is  obtained  if  the  non- trivial  solution 
is  found  from  equation  (92).  It  is  to  be  noted  that  the  buckling  leads  in 
fases  1 to  5 can  all  be  obtained  from  eq.  (107),  as  they  should  be. 

If  we  define  the  stress  ratios  R^, 
jng  formulas 


Rg,  R^,  according  to  the  f'oilov- 


\ - — » 

C Buckling  stress  under  compression  alone 


Critical  compressive  stress 


Critical  bending  moment 


IC  Buckling  moment  under  bending  alone 


2fra 


Critical  torsional  moment 

Buckling  moment  under  torsion  alone 


then  equation  (105)  may  be  written  as 


R>P  " - \/l  - (Rg  + Rc). 


(108) 


The  interrelationship  between  compression,  bending,  and  torsion- 
al stress  ratios  given  by  equation  (108)  is  plotted  for  engineering  use  in 
Figs,  10  and  11.  Once  any  two  stress  ratios  are  specified,  the  buckling 
value  of  the  remaining  stress  ratio  can  be  determined  graphically  from 

fV.oQO  enwa  a 
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Experimental  Inv eg  Ligation* 
Case  1.  Buckling  Under  Axial  Compression 


Test  program. 

In  an  effort  to  test  cylinders  within  a range  of  practical 
interest,  it  is  desirable  that  the  radius  to  thickness  ratio  should  be  as 
large  as  possible.  Whereas  a relatively  large  value  of  the  a/t  parameter 
is  easily  obtainable  for  a homogeneous  cylinder  by  constructing  the  cylinder 
of  sheet  of  several  thousandths  of  an  inch  thickness,  the  sandwich  cylinder 
presents  the  problem  that  considerably  greater  thickness  is  required  for 
the  sandwich  for  practical  reasons, 

A l/B-in,  thick  sandwich  was  considered  the  minimum  thickness 
feasible,  and  consequently  rather  large  radii  were  chosen  so  that  the  Jt/t 
values  would  be  within  a range  of  possible  use  in  aircraft  construction. 
Three  different  radii  were  used  to  permit  a variation  in  the  0/t  parameter 
to  be  studied  within  the  capacity  of  the  testing  facilities.  To  determine 
the  length  effect,  if  any,  two  different  lengths  were  used  in  the  investiga- 
tions. 


All  other  details  of  the  specimens  were  similar.  The  24S-T  alu- 
minum-alloy faces  of  0. 010-in,  thickness  were  bonded  to  a l/8-in.  cellular 
cellulose  acetate  (CCA)  core  of  approximately  4.5  to  5*0  lbs,  per  cu.  ft. 
density.  At  each  end  of  the  cylinder,  a 1 l/2-in,  deep  hardwood  insert 
replaced  the  core  to  allow  for  the  uniform  introduction  of  end  loads  into 
the  cylinder* 


The  cylinders  were  manufactured  by  Skyuyne,  Inc,,  Port  Jervis, 

N.Y.,  and  their  fabrication  details  are  quoted:  "All  aluminum  faces  are  coated 
with  a metal  primer  prior  to  assembly  to  the  core.  The  primer  is  an  adhesive 
manufactured  by  the  Casein  Company  of  America  and  bears  their  code  number 
of  NT-442. 


* The  experimental  work  was  carried  out  mostly  by  G*  Gerard  and  F0  K. 
Teichnann. 
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The  core  is  then  planed  to  the  thickness  of  l/S  in0,  and  solid  wood  re- 
inforcements are  glued  to  the  edge.  A secondary  adhesive  is  used  to  combine 
the  aluminum  fac^s  and  the  corej  the  adhesive  used  is  Plastron  250-2.  Al- 
though this  resin  adhesive  will  polymerize  at  a room  temperature  of  75°?oV 
sufficient  heat  was  applied  during  the  fabrication  of  the  cylinders  to  ai- 
cslerate  the  "setting”  of  the  bond.  Pressure  necessary  to  bring  the  surface.'- 
in  intimate  contact  was  applied,  and  the  assembly  was  molded  in  a period 
of  2 to  3 hours.  The  final  operation  consists  of  trimming  or  finishing  the 
ends  of  the  cylinder  to  assure  an  even  plane  and  squareness  to  the  ends0rt 

Load  was  applied  to  each  spociraent  by  means  of  a Baldvin-Southwark 
universal  hydraulic  testing  machine.  The  10-in.  radius  cylinders  were  tested 
ir.  a machine  of  60, 000-lb.  capacity,  whereas  the  remaindor  were  tested  in 
a 200,000-lb.  capacity  machine.  The  specimen  rested  on  a l/4-in.  machined 
plate,  which  in  turn  rested  directly  on  the  base  of  the  machine.  A 1 1/2- 
in.  machined  steel  plate  was  centered  on  top  of  the  cylinder  and  load  was 
applied  to  this  through  a self-aligning  compression  head. 

During  the  progress  of  the  test,  load  was  applied  slowly.  At 
appropriate  increments  of  load,  strain-gage  readings  at  four  stations  and 
deflection  re  ad  in  gp  at  three  stations  were  taken.  Near  the  buckling  load 
of  the  cylinder,  the  load  increments  were  reduced.  At  failure,  both  the 
maximum  load  and  the  drop  load  at  which  the  machine  stabilized  were  recorded. 

ultimate  loads  and  dimensions 

A summary  of  the  important  test  results  is  given  in  Table  2. 
The  various  dimensions  given  in  this  table  are  those  that  were  obtained 
prior  to  test  at  the  estimated  center  of  the  buckle  pattern.  It  was  noted, 
however,  that  this  location  usually  did  not  coincide  with  the  location  of 
minimum  thickness  or  maximum  length  of  the  cylinder. 


Table  2 


Presentation  of  Teat  Data 


Cylinder* 

Nominal 

Length, 

In. 

Inside 
Rac  .us, 
In# 

Total 

Thickness, 

In. 

Face 

Thickness, 

In, 

Calculated 

Core 

Thickness, 
In,  . _ 

Ultimate 

Load, 

Lbs. 

Drv>p 

Load, 

Los*. 

1C12A+ 

12 

lo.o*; 

0.144 

0.0107 

0.1226 

9,675 

4,22.? 

1012B 

12 

9.99 

0.134 

0.01(77 

0.1126 

31,825 

4nC2f 

1024A 

24 

10.09 

0,138 

0.0107 

0.1166 

10,625 

5*025 

1024B 

24 

10.05 

0.135 

0.0107 

0.1136 

13,985 

5,625 

1212A 

12 

12.00 

0.142 

0.01(77 

0.1206 

18,175 

5,675 

1212B 

12 

11.89 

0.136 

0.0107 

- 0.1146 

18,475 

3.950 

1224A 

24 

11.88 

0.132 

0.01C7 

0.1106 

17,625 

8,225 

1224B 

24 

Premature  failure  due  to 

poor  bonding 

1412A 

12 

13*89 

0.143 

0.0107 

0.1216 

26,950 

2,250 

1412B 

12 

13.86 

0.134 

0.0107 

0.1126 

33,800 

2,400 

1424A 

24 

13.88 

0.132 

0.01C7 

0.1106 

21,150 

5,900 

1 / O/  T> 

24 

13.84 

0.131 

0.0107 

0.1096 

18,450 

1 £ 
8 

* Note  that  the  cylinders  are  identified  by  the  radius  that  is  given  by  the 
first  two  digits  and  the  length  that  is  given  by  the  last  two  digits# 


+ This  cylinder  was  initially  buckled,  dicing  manufacture,  at  the  location  cf 
buckling  under  load. 


Material  property  tests  conducted  on  the  CCA  core  material  in- 
dicated that  a shear  modulus  of  Gq  ■ 2,000  lbs.  per  sq*  in*  - 20  per  cent 
was  a representative  value*  Since  the  stress  in  the  faces  at  buckling 
was  well  below  the  proportional  limit  of  the  material,  no  material  property 
tests  were  conducted  on  the  face  material# 

Description  of  failure 

In  general,  the  mode  of  instability  observed  on  the  cylinders 
tested  can  be  roughly  divided  into  three  groups J (a)  Cylinders  in  which 
buckles  appeared  in  isolated  locations  before  snap  buckling  precipitated 
failure;  1012A,  1012B,  1212B.  (b)  Cylinders  in  which  no  visible  evidence 


of  buckling  occurred  before  snap  buckling  precipitated  failure:  1024A, 

1024B,  1212A,  1224B.  (c)  Cylinders  in  which  outward  bulging  either  pre- 

ceded Or  accompanied  snap  buckling:  1412A,  1412B,  1424A9  L424B. 

The  cylinders  in  group  (a)  exhibited  isolated  buckles  that  grew 
as  the  load  was  increased.  These  buckles  usually  formed  the  nucleus  of 
the  snap  buckle  pattern  that  occurred  at  failure  with  explosive  violence. 

At  failure,  the  ultimate  load  dropped  to  a value  approximately  40  per  cent 
of  the  peak  value.  The  behavior  of  group  (b)  was  simil  =r  with  respect 
to  the  drop  in  load.  The  characteristic  buckle  pattern  obtained  in  these 
groups  is  shown  in  Fig.  12. 

Group  (c)  cylinders  exhibited  a short  wave  length  outward  bulge, 
usually  at  the  junction  of  the  CCA  core  and  hardwood  insert,,  in  addition 
to  a snap  buckling  pattern.  The  former  mode  of  instability  appeared  to 
be  characteristic  of  the  type  associated  with  wrinkling  of  the  faces.  In 
the  cases  where  the  otttward  bulging  preceded  failure,  it  was  noticed  that 
portions  of  the  bulge  area  changed  to  the  character-*  **tic  snap  buckle  pattern 
at,  or  immendiately  after,  failure.  In  addition,  the  drop  load  was  much 
less  than  for  the  other  two  groups,  averaging  about  20  per  cent  of  the 
peak  value. 


Buckle  patterns 

Data  on  the  buckle  patterns  obtained  in  the  tests  ars  giver: 
for  °ach  cylinder  in  Table  3. 

The  shapes  of  the  buckles  obtained  from  the  .snap  buckle  patterns 
were  diamond-shaped  as  shown  in  Fig.  12  and  of  approximately  3 to  1 aspect 
ratio  as  given  in  Table  3.  They  indicated  the  characteristic  tendency 
to  buckle  inwardly.  It  is  noted  for  purposes  of  comparison  that,  for  a 
homogeneous  cylinder,  the  buckle  is  usually  diamond-shaped  and  of  an  aspe,  t 
ratio  close  to  unity. 

Another  detail  r'  importance  was  the  occurrence  of  only  one 
buckle  in  the  axial  direction  for  all  the  cylinders.  In  testing  Cylinders 
1012B  and  1024A,  excessive  deformation  was  applied  to  the  cylinder  tc  d«  ■ 


T 
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termine  if  additional  buckles  would  form  in  the  axial  direction.  It  was 
observed  that  no  further  buckling  occurred.  'This  behavior  is  probably 
due  to  the  formation  of  a weak  area  at  the  buckle  location  once  buckling 
occurs.  Further  deformation  of  the  cylinder  appears  to  be  completely 
restricted  to  the  buckle  location  with  ultimate  crushing  of  the  cere  and 
separation  between  the  face  and  core. 

Table  3 

Data  on  Buckle  Patterns 
Buckle  Size 

Longitu=  Circum-  Number  of 


Cylinder 

dinally 
= i/m 
in. 

ferentially 

LO  = ZS 
n 

in. 

Circum- 

ferential 

Buckles* 

LO  X. 

T 

"x  n 

1012A 

1.5 

4.0 

24 

2.47 

1P12B 

1.5 

4.0 

24 

2.67 

1024A 

1.5 

4.0 

24 

2.67 

1212A 

1.5 

4.5 

24 

3.00 

1224A 

1.5 

4.5 

24 

3.00 

* Note  that  these  buckles  were  sinusoidally  tu- ranged  around  the  circum- 
ference (see  Fig.  1|)» 

Influence  of  length 


A 12-and  24-in » length  of  each  cylinder  were  included  in  the 


test  program  tci  determine  trie  influence  of  .Length  upon  the  buckling  char- 
acteristics within  this  range  of  length.  Allowing  for  the  hardwood  Inserts 
in  the  end  of  each  specimen,  the  free  length,  of  the  cylinders  in 
which  buckling  is  likely  to  occur  Is  reduced  to  9 ana  21  in.,  respectively. 


The  buckling  load  per  inch,  N ^ as  a function  of  Lr/%  is  shown, 
in  Fig.  1JL  Within  this  range  of  the  parameter  and  considering  the  sCuttex 
of  the  test  data,  no  conclusions  as  to  the  effect  of  length  could  be  de- 
termined. Examination  of  the  buckle  dimension  data  indicates,  however, 
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that  no  affect  of  length  on  the  buckling  load  should  occur,  since,  for 
the  shortest  cylinder , the  axial  wave  leagth  was  approximataly  onc-sia-th 
of  the  free  length.  It  is  anticipated  that  only  when  the  free  length 
approaches  the  ucial  length  will  the  length  influence  the  buckling  char- 
acteristics. 

Physical  and  material  property  variation 

Examination  of  the  dimensional  properties  of  each  cylinder 
indicated  that  the  variation  in  dimensions  from  the  average  were  of  the 
following  approximate  orders  of  magnitude:  (a)  face  thickness,  t - 0.0015 
in.  (commercial  tolerance)}  (b)  core  thickness,  h - 0.008  in.;  (c)  inside 
radius;  a - 0.030  in.;  (d)  length,  1 - 0.008  in. 

Variation  in  the  shear  modulus  of  the  CCA  core  was  found  to 

be  G_  - 20  per  cent.  The  variation  of  E..  of  the  face  material  is  known 

c 1 6 

to  negligible  for  aluminum  alloys.  A value  of  Ef  * 10.5  x 10  lbe. 

per  8q0  in.  was  used. 

An  attempt  was  made  to  correlate  the  location  estimated  to  be 
responsible  for  buckling  with  initial  imperfections  such  as  minimum  values 
of  thickness  or  maximum  values  of  length.  From  the  tabulation  of  dimen- 
sional data  for  each  of  the  cylinders,  no  correlation  was  indicated. 

Since  correlation  with  dimensional  imperfections  was  lacking,  it  is  assum- 
ed that  the  buckle  location  is  caused  primarily  by  local  weak  spots  created 
by  the  variations  in  the  core  shear  modulus.  The  fact  that  only  one  Longi- 
tudinal buckle  formed  in  all  cylinders  appears  to  be  further  conf  rmation 
of  this  hypothesise 

Additional  evidence  of  this  type  of  behavior  of  sandwich  elements 

1? 

subject  to  buckling  has  previously  been  offered  by  Hoff  and  Hammer 
They  suggest  that,  in  a sandwich  structure  subject  to  uniform  load,  such 
as  a beam  in  bending,  the  results  are  a statistical  average  of  the  core 
pi  ujijbi'tloa.  On  tlio  ulLor  hand 9 a a true  Lure  subject  to  Duelling  responds 
to  local  veak  spots  caused  by  defects  in  the  core  properties,  and,  hence, 
the  results  should  be  based  on  the  lowest  properties  rather  than  cm  the 
average  properties. 
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Correlation  of  theory  and  test  data 

It  was  found  herein  that  for  values  of  E^t/G^*  ^ 0.95»  the 
critical  loading  is  given  by  the  following  equation? 

Ncr  - (h  + t)  Gc 

It  appears  that,  for  this  range  of  E^t/GQV  the  critical  stress  is  inde- 
pendent of  the  radius  l^and  depends  essentially  upon  the  core  properties. 

To  test  the  validity  of  this  conclusion,  the  test  data  are  plotted 
in  Fig.  as  a function  of  ft.  Also  shown  are  the  results  obtained  from 
the  preceding  equation  based  on  the  upper  and  lower  limits  of  the  shear 
modulus,  0 , given  previously.  Since  no  correlation  of  failure  location 
with  dimensional  imperfections  was  obtained  as  discussed  in  the  previous 
section,  it  is  assumed  that  the  variation  in  the  shear  modulus  is  primarily 
responsible  for  buckling. 

While  the  test  data  in  Fig.  14  appear  to  indicate  an  apparent 
influence  of  the  radius,  it  is  **elt  that  this  apparent  effect  can  be  ex- 
plained by  the  variation  in  thf  shear  modulus.  Consideration  must  also 
be  given  to  the  fact  that  the  10-in.  radius  cylinders  were  tested  in  a 
60, OOO-lb.  capacity  testing  machine,  whereas  the  remainder  were  tested 
in  a 200,000-lb.  machine.  TMs  may  have  some  bearing  on  the  lower  values 
of  the  10-in.  radius  group.  Furthermore,  the  cylinders  were  manufactured 
in  order  of  increasing  radius,  and,  if  some  allowance  is  made  for  an  im- 
provement in  manufacturing  technique  with  experience,  then  the  large  radius 
cylinders  are  favored. 

Considering  all  these  factors,  it  appears  that  reasonable  corre- 
lation between  the  critical  load  obtained  from  the  theory  and  the  test  data 
was  obtained  within  the  region  of  E^t/Gc a of  the  experiments. 


Case  20  Buckling  under  torsion 


Test  program 

Again  the  sandwich  used  had  a thicknees  of  i/8  inch0  The 
nominal  radius  of  6 inches  which  was  used  for  all  cylinders  was  dictated 
by  the  capacity  of  the  testing  facilities#  To  determine  the  effect  of 
length  upon  the  buckling  load,  two  different  lengths  were  used  in  the  in- 
vestigation# 

The  most  significant  parameter  in  overa.';.!  buckring  of  sandwich 
elements  is  the  shear  rigidity  of  the  cere  material#  Consequently,  two 
different  materials  were  selected  for  the  cores  cellular  cellulose  acetate 
(CCA)  of  approximately  5#  5 to  6 lbs/cu#  ft#  density  and  end  grain  balsa 
(EG®)  of  approximately  5 to  5®5  lbs/cu.  ft#  density#  For  all  cylinders, 
the  face  material  used  was  0o010  inch  24S-T3  aluminum  alloy# 

To  faciltate  instalifcion  of  the  test  specimen  in  the  loading 
jig,  two  inch  wide  hardwood  blocks  were  bended  to  both  the  inside  and  out- 
side faces  at  each  end  of  the  cylinder  (see  Fig#  15)#  The  end  faces  of  these 
blocks  were  then  turned  down  on  a lathe  to  insure  planeness  of  the  face, 
perpendicularity  with  the  axi3  of  the  cylinder  and  parallelism  of  the  faces# 

Test  procedure 

The  pertinent  dimensions  of  each  specimen  were  determined 
as  follows.  Four  evenly  spaced  longitudinal,  stations  were  lightly  scribed 
on  the  outside  wall  of  tbs  cylinder#  The  outside  radius  was  then  measured 
to  tha  nearest  0.01  inch  at  each  pair  of  stations  at  both  ends#  The  length 
of  the  cylinder  between  hardwood  blocks  was  also  determined  at  each  station 
to  the  nearest  G#01  inch# 

The  total  thickness  of  the  cylinder  wall  was  measured  at  each 
station  at  both  ends  just  inside  the  hardwood  blocks  and  at  the  center  by 
means  of  a special  dial  gage  with  a least  division  of  OoOOl  inch.  The 
thickness  of  the  ft|c«s  was  measured  with  a lengthens and ths  micrometer  before 
the  cylinders  were  assembled.  The  average  dimensions  for  each  cylinder  are 
given  in  Table  4* 
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Installation  of  the  specimens  in  the  specially  designed  torsion 
loading  jig  was  accomplished  by  bolting  the  specimen  to  the  steel  end 
plates  through  the  hardwood  blocks  at  the  ends  of  the  cylinder.  The  loading 
mechanism  is  shewn  in  Fig.  15  with  a specimen  installed  and  tested,.  Load 
was  applied  by  means  of  a hydraulic  tension  jack  and  was  measured  by  a 
calibrated  strain  gage  tension  link  at  the  jack  which  was  accurate  to  - 
2 per  cent. 


Physical  property  tests 

The  important  physical  properties  in  overall  torsional  buck- 
ling of  sandwich  construction  are  the  shear  stress-strain  characteristics 
of  the  face  material  and  the  shear  modulus  of  the  core  material.  The  methods 
used  to  determine  these  properties  are  discussed  in  some  detail  due  to 
their  importance  in  evaluation  of  the  teat  data. 

The  shear  moduli  of  the  core  materials  were  determined  by  test- 
ing as  simple  beams,  strips  of  flat  sandwich  construction.  The  beams  were 
fabricated  by  the  manufacturer  of  the  cylinders  of  the  same  materials  and 
in  the  3ame  manner  as  the  cylinders  themselves.  The  test  consisted  of  ob- 
taining load-deflection  data  from  which  the  shear  modulus  of  the  core  mater- 
ial was  determined  from  the  difference  in  deflection  between  the  experi- 
mental valus  and.  that  computed  for  the  beam  neglecting  shear  deflections. 

From  six  tests  of  each  core  material,  the  following  values  of  shear  modu- 
lus were  obtained: 

cellular  cellulose  acetate:  G « 3400  psi  - 20$ 

c 

end  grain  balsa:  Gn  = 15,700  psi  - 7$. 

No  experimental  technique  is  known  for  d irectly  determining  shear 
stress-strain  characteristics  of  thin  sheet  material.  The  method  used  in 
this  investigation  was  to  obtain  the  stress-strain  properties  under  axial 
comoression  and  tension  loads  of  specimens  oriented  at  45  degrees  with  the 
grain  direction  of  the  face  material.  This  orientation  was  chosen  since 
a put's  shear  field  can  be  represented  by  orthogonal  tension  and  compression 
stress  of  the  same  magnitude  as  the  shear' stress  at  45  degrees  with  the 
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dircction  of  the  shear  stress.  The  axial  stre3s-strain  data  so  obtained 
are  shown  in  Fig,  16  aa  the  limits  of  the  variation  in  stress-strain  prop- 
erties from  the  average  of  the  specimens  tested, 

A shear  stress-strain  curve  can  be  constructed  from  axial  stresa- 
strain  data  by  use  of  the  following  transformations  of  the  maximum  shear 
theory* 

/ 

t “ o/2  and  X *=  (l  + v)e  Q09) 

In  the  yield  region  of  the  stress-strain  curve,  Poisson's  ratio  can  be 

. , ,13 

v » 0.5  - 0.2  E^/E  (110) 

where  E and  E are  the  elastic  and  secant  moduli  of  the  axial  strasc- 
a 

strain  r_  <m-  The  desired  shear  stress-strain  curve  shown  in  Fig.  l6  was 
obtained  by  applying  Eqe.  (109)  and  (110)  to  the  upper  and  lower  limits 
of  the  axial  stress-strain  data. 

Test  results 

Average  measured  dimensions  of  the  sandwich  cylinders  to- 
gether with  the  failing  loads  are  given  in  Table  4.  Photographs  of  typical 
buckle  patterns  t.eVon  after  conclusion  of  the  tests  are  shown  for  each 
type  of  cylinder.  Fig.  17  shews  an  end  grain  balsa  core  cylinder.  Fig. 

18  shows  a CCA  cylinder  of  24  inch  length  and  Fig.  19  a CCA  cylinder  of 
12  inch  length. 

It  is  noted  that  in  all  cases,  there  was  no  observation  of  buck- 
ling before  the  failing  load  was  reached. 

Correlation  of  theory  with  experimental  results 
(a)  End  grain  balsa  cylinders 

For  this  group  of  cylinders,  we  find  that  the  parameters 
E-pt/G^a  = 1.1  and  a/(h+t)  * 48.  These  cylinders  behave  essentially  as 
homogeneous  cylinders.  Therefore,  the  results  of  Ref.  14  can  be  used  to 
evaluate  the  effect  of  the  boundary  conditions  and  leng  li  of  the  cylinder 
upon  the  critical  shear  stress.  From  Ref.  14,  for  a cylinder  clamped  at  the 
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ends 


t (2t)l2 
cr 

tt2!) 


n q^z 

WV  y 


3/4 


(111) 


where  D ■ bending  rigidity  of  sandwich  plate  ■ >‘^t ( h+t ) ~/2 ( l-v^) 

2 1/2  2 

Z “ curvature  paremeter  = (l-v^,)  i-fc/2at 

In  using  the  above  relations  for  e sandwich  plate,  the  thickness  of  the 
homogeneous  cylinder  has  been  replaced  by  twice  the  face  thickness  of  the 
sandwich  cylinder,  (2t);  Eq.  (ill)  can  be  reduced  to  the  following  by  taking 
v = 0.3. 

i 3/ A 

Tcr  - 1.44Ef(h+t)2A^"  (at)>  (112) 


By  substituting  the  appropriaxe  values  given  in  Table  4 for  the  parameters 
appearing  in  Eq.  (112)  i+  is  found  that 

t - 433,000  psi  (113) 

cr 

T^e  value  of  the  critical  stress  given  by  Eq.  (113)  is  based  on  elastic 
buckling  and  consequently  is  far  in  excess  of  the  shear  proportional  limit 
of  the  material  ijsed.  A method  for  computing  plastic  shear  buckling  stress- 
es is  given  in  Ref*  15,  in  which  it  is  suggested  that  the  effective  modulus 
to  be  used  in  plastic  shear  buckling  is  the  shear  modulus.  By  this  method, 
the  critical  shear  strain  is  computed  and  the  plastic  buckling  stress  is  then 
obtained  from  the  shear  stress-strain  curve  at  the  particular  value  of 
critical  strain.  From  Eq.  (113),  therefore, 

r 

Y = T /„  » 0.109  (11/,) 

‘'cr  cr/G^, 

This  is  evidently  a very  large  value  of  critical  sorain  for  which 
the  shear  stress-3train  cvr^ve  is  not  known.  In  addition,  it  is  highly 
doubtful  that  the  method  given  in  Ref.  15  (which  is  for  the  yield  region 
of  the  stress-strain  curve)  would  apply  at  such  large  values  of  critical 
strain.  Therefore,  to  effect  a comparison  between  theory  and  test  data, 
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the  test  date  are  arbitrarily  plotted  in  Fig.  21  at  a strain  of  0.01  to  ob- 
serve "che  relation  between  these  date  and  the  shear  stress-strain  our/e 
given  in  Fig.  16.  These  data  are  plotted  at  s strain  of  0.01  since  this 
appears  to  be  a reasonable  value  indicating  the  limit  of  the  yield  region 
of  th9  stress-strain  curve.  Although  reasonable  agreement  is  indicated 
in  Fig,  21  it  is  to  be  noted  that  significant  difference  between  the  theory 
and  tests  may  be  obscured  by  the  ffact  that  buckling  occurred  in  the  plastic 
r0gion«  , 


Tabl  e 4 

AVERAGE  DIMENSIONS  AND  FAILING  LOADS  OF  SANDWICH  CYLINDERS 


Cylin- 

der 


Core 

Materi- 


Total  Average 

Ffean  Thick-  Face  Failing 

TH  oje— . 

Radius  ness  ness  Moment 


Failing  Failing 
Load/in.  Stress 


1 

a 

(h+t) 

t 

Mcr 

Ncr 

^cr 

No.l 

EGB 

20.06" 

5.99" 

0.145" 

0.103" 

l26,OC0"lb, 

560  lb/in. 

28,000psi 

No, 2 

EGB 

19.50 

5.99 

0.155 

0.103 

103,600 

460 

22,900 

No.  3 

CCA 

7.90 

6.01 

0.138 

0.105 

42,000 

185 

8,550 

No.  4 

CCA 

7,95 

6.00 

0.142 

0.105 

25,200 

111 

5,320 

No.  5 

CCA 

19.80 

6.01 

0.137 

0.105 

26,400 

117 

5,550 

No.6 

CCA 

19.92 

6.02 

0.132 

0,105 

22,700 

100 

4,730 

(b)  Cellular  cellulose  acetate  cylinders 


The  average  value  of  the  parameter  Eft/G  a for  this  group  of  cylin- 
ders  was  approximately  five.  According  to  Fig.  20,  the  buckling  strength 
of  tills  group  of  cylinders  should  be  effected  only  by  the  shear  rigidity  of  the 
core?  Therefore,  the  experimental  values  of  the  torsional  instability 
coefficients  were  computed  as  shown  it;  Table  5 and  plotted  in  Fig.  20. 


I 
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Table  5 

EXPERIMENTAL  TORSION  INSTABILITY  COEFFICIENTS 


FOR  CCA  CYLINDERS 

fi-t/G  a.  = 

5,'i  5 J 

^\h+t) 

a 48 

Cylinder 

N 

cr 

h+t 

°e 

K.-N  /{h 
t.  cr 

No.  3 

185  lb. /in. 

0.118" 

3,400  psi 

0.462 

No.  4 

Ill 

0.122 

3,400 

0.268 

No.  5 

117 

0.117 

3,400 

0.294 

No.  6 

100 

0.112 

3,400 

0.263 

The  comparison  of  the  theory  and  test  data  shown  in  Fig.  21  re- 
veals that  the  test  results  are  considerably  below  the  theoretical  vi  lue. 
The  large  variation  in  shear  modulus  of  the  core  (-  20%)  is  certainly  not 
sufficient  to  account  for  this  discrepancy.  It  appears,  however,  that 
poor  bonding  of  the  faces  to  the  core  may  be  the  cause  of  the  discrepancy. 

It  is  possible  that,  because  of  the  curvature  of  the  structure, 
inspection  by  tapping  the  cylinders  did  not  reveal  any  area  of  poor  banding 
before  the  test.  During  loading,  however,  evidences  of  poor  bonding  were 
noticed.  The  evidence  was  a crackling  noise  which  was  heard  almost  immedi- 
ately upon  the  application  of  load. 


Additional  evidence  was  obtained  after  the  cylinders  were  tested. 
In  some  of  the  cylinders,  the  faces  were  peeled  from  the  core  and  it  was 
■c'ound  that  there  were  considerable  areas  in  which  no  core  material  remained 
attached  to  the  bonding  of  the  faces.  It  appears  that  in  a satisfactory 
bond,  particles  of  the  core  materia],  would  adnere  bo  the  bonding  of  the 
faces  after  peeling. 


This  difficult;'-  occurred  only  with  the  CCA  core  material  and 
may  possibly  be  attributed  to  the  small  thickness  (l/8")  of  core  material 
used.  In  addition  the  CCA  was  produced  by  a new  manufacturer  and  there 
was  evidence  of  considerable  nonuniformity  of  cell  size. 
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Case  3 r Buckling  under  bending 
Test  program 

The  cylinders  tested  in  pure  bending  had  again  a thickness 
of  1/8  inch  and  a radius  of  6 inches.  To  ensure  uniformity  of  introduc- 
tion of  bending  load,  all  the  cylinders  were  36  inches  long.  Two  different 
material  were  again  used  for  the  coral  cellular  cellulose  acetate  (CCA) 
of  approximately  6 lbs./cu.  ft.  density  and  end-grain  balsa  (EGB)  of  approx- 
imately 5 to  5»5  lbs./cu.  ft.  density.  For  all  cylinders,  the  face  material 
used  was  0.010  inch  24S-T3  aluminum  alloy. 

To  faciltate  installtion  of  the  test  specimen  in  the  loading 
jig,  two  inch  wide  hardwood  blocks  were  bonded  to  both  the  inside  and  out- 
side faces  at  each  end  of  the  cylinder  (see  Figs.  22  and  2Z.).  The  end  faces 
of  these  blocks  were  then  turned  down  on  a lathe  to  insure  planeneas  of  the 
face,  perpendicularly  with  the  axis  of  the  cylinder  and  parallelism  of  the 
faces. 

Test  procedure 

The  pertinent  dimensions  of  each  specimen  were  determined  as  in 
the  case  of  buckling  under  torsion.  The  average  dimensions  for  each  cylin- 
der are  given  in  Table  6. 

Installation  of  the  specimens  in  the  specially  designed  banding 
loading  jig  vas  accomplished  by  bolting  the  specimen  to  the  steel  end  plates 
through  the  hardwood  blocks  at  the  ends  of  the  cylinder.  The  loading  jig 
is  shown  In  Fig.  op  with  a specimen  installed  and  tested.  Load  was  applied 
by  metns  of  a hydraulic  tension  jack  and  was  measured  by  a calibrated  strain 
gage  tension  link  at  the  jack  which  was  accurate  to  - 2 percent. 

Physical  property  tests 

The  important  physical  properties  in  overall  bending  buckling 
of  sandwich  construction  are  the  compressive  stress-strain  characteristics 
of  the  face  material  and  the  shear  modulus  of  the  core  material. 
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The  with-grain  compressive  stress-strain  characteristics  of  the 
0.010  in.  24S-T3  aluminum  alloy  used  for  the  faces  of  the  sandwich  cylinders 
was  determined  by  use  of  a sol -Id  guide  compression  jig  of  the  National. 
Bureau  ■ • ■ • a.  . ,*\  pu . ino  2tjX"eE?“— surair*  oax>a  are  shown  in  Fxg.  23. 

The  shear  modjli  of  the  core  mace rials  were  determined  by  testing  as  be- 
fore. From  at  least  four  tests  of  each  core  material,  the  following  values 
of  shear  modulus  were  obtained! 


cellular  cellulose  acc+«tej  G„  * 3400  pci  - 17? 
end  grain  balsa:  G_  m 17,600  psi  - 15? 


Average  measured  dimensions  of  the  sandwich  cylinders  together 
with  the  failing  (buckling)  loads  are  given  in  Table  6.  Fig.  24-  shows  a 
typical  overall  buckle  pattern  obtained  on  most  of  the  1/8-in.  cylinders j 
additional  data  on  the  buckling  patterns  are  given  in  Table  7.  The  measured 
wavelengths  were  obtained  in  regions  of  relatively  little  distortion. 

Table  6 

AVERAGE  DIMENSIONS  AND  FAILING  LOADS  OF  SANDWICH  CYLINDERS 


Cylin- 

der 

Core 

Matsri 

Total 

Mcsri  Thick- 
al  Length  Radius  ness 

Average 

Face 

Thick- 

ness 

Failing 

Moment 

Maxi- 

mum 

Failing 

Load/iru 

Maxi- 

mum 

Failing 
St  re  35 

i 

a 

(h+2t) 

t 

ncr 

N mH/v&~ 

or 

0 

cr 

No.l 

CCA 

32.09" 

6.48" 

o.ui» 

0.0105" 

67 ,900»lbs . 5141t/in. 

19,300  psi 

No  z? 

CCA 

32.13 

4 rn 

JKJ 

0-125 

0.0105 

46,100 

347 

16,500 

No. 3 

CCA 

32.13 

6.49 

0,125 

0.0105 

77,700 

587 

28,000 

No.4 

CCA 

32.16 

6.53 

0.141 

0.0105 

89,600 

670 

31,900 

No.  5 

EGB 

31.29 

6.47 

0.141 

0.0105 

53,100 

404 

19,300 

No.6 

EGB 

31.09 

6.55 

0.156 

0.0105 

65,700 

500 

23,800 
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7 able  7 


DATA 

ON  INSTABILITY  PATTERNS 

TVpe  of 

Longitudinal 

Circumferential 

Cylinder 

Instability 

Location 

Wave  length 

Wave  Length 

L 

X 

L© 

No.l  l/8"CCA 

Buckle 

Central 

1" 

3" 

No. 2 1/8  "CCA 

Buckle 

Central 

1 

** 

No. 3 1/8  "CCA 

Buckle 

Centr  al 

1 1/4 

3 1/2 

Antisymmetric al 

No. 4 1/8  "CCA 

Wrinkle 

Both  end  s 

7/16 

ira 

No. 5 l/8»EGB 

Buckle 

Central 

1 1/4 

3 1/2 

Antisymmetrical 

No. 6 1/8 "EGB 

Wrinkle 

At  end 

7/16 

TT& 

4 

Oerrel-rfcion  of  Theory  with  Hacpwrinental  Results 

The  l/8  in.  CCA  cylinders  are  sandwich  cylinders  with  weak 
cores.  The  average  centroidal  height,  (h+t),  of  cylinders  No.  1 to  No.  4 
is  0.118  in.  - 6%  and  the  average  shear  modulus  is  3400  psi  - 17%.  There- 
fore, the  average  shear  rigidity  of  this  group  of  cylinders  is  400  lbs./in. 
+ 23%,  which  is  the  theoretical  buckling  loading.  The  average  experimental 
failing  loading  was  found  to  be  530  lbs./in.  - 33%.  Thus,  the  experimental 
failing  loading  is  32  percent  higher  than  the  theoretical  value,  on  the 
average. 

The  fact  that  the  bending  buckling  load  is  higher  than  the 
compressive  buckling  load  for  cylinders  of  the  same  din»nsions  has  been 
observed  in  tests  on  homogeneous  cylinders.  It  has  been  suggested  that 
this  behavior  may  be  due  to  the'  buckling  phenomena  responding  to  the  aver- 
age value  of  compressive  stress  on  the  cross  section  rather  than  the  maxi- 
mum stress.  In  the  case  of  pur  bending,  the  maximum  stress  would  be  1.4 
t-imes  £3  great  as  the  average  value. 

For  sandwich  cylinders  in  compression,  the  theory  indicates 
tiiac  tne  value  oi  the  parameter  E^t/G^a  governs  the  behavior  of  the  cylin- 
der. Thus,  it  has  been  found  that  cylinders  can  be  considered  weak  in 
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shear  when  the  value  of  E^t/Gca  exceeds  unity.  To  effect  a comparison 
between  the  cylinders  loaded  in  compression  and  bending  which  were  differ- 
ently constructed  as  concerns  Gc  and  a,  it  appears  that  if  the  values  of 
E^t/G^a  were  equal  or  at  least  greater  than  unity  for  each  group,  a direct- 
comparison  would  be  valid. 

The  value  of  this  parameter  for  the  compression  specimens 
averaged  5 and  for  the  bending  specimens  E^t/Gca  ■ 5 also.  Thus,  since 
the  compressive  test  data  checks  theory  very  veil,  it  can  be  concluded 
that  in  bending  the  32?  higher  experimental  value  is  In  agreement  with  the 
behavior  observed  for  homogeneous  cylinder*-  thus  accounting  for  the  dis- 
crepancy. 

Although  these  conclusions  are  based  on  data  with  which  large 
percentage  variations  are  associated,  it  appears  that  the  trend  is  clearly 
indicated.  The  fact  that  the  scatter  of  the  experimental  valuos  is  con- 
siderably higher  than  that  for  the  core  properties  is  based  on  the  fact  that 
the  latter  values  were  established  from  beam  bending  tests.  In  a beam 
test,  the  results  are  a statistical  average  of  the  core  properties,  whereas 
a structure  subject  to  buckling  responds  to  local  weak  spots  caused  by 
defects  in  the  core  properties.  Thus,  the  scatter  would  tend  to  be  greater 
for  the  buckling  data. 

The  l/8  in,  end  grain  balsa  cylinders.  No.  5 and  No.  6,  have 
cores  much  stronger  in  shear  than  the  cylinders  with  CCA  cores.  Based  on 
the  observations  made  on  these  two  cylinders,  it  appeared  that  they  cannot 
be  considered  in  ths  weak  in  shear  category.  These  cylinders,  therefore, 
failed  at  a small  fraction  of  the  ■theoretical  load.  This  is  evidenced  by 
the  fact  that  the  average  falling  loading  of  this  group  was  452  Ibs./in* 
as  compared  to  5"0  lbs./in«  for  the  1/6  in. GCA  cylinders,  although  the  shear 
modulus  of  the  EGB  core  was  over  five  times  as  great  as  the  GCA.  core.  The 
fact  that  the  EGB  cylinders  were  not  even  as  strong  as  the  CCA  cylinders  is 
difficult  to  reconcile  with  the  core  properties.  Unless  some  other  mode 
of  buckling  is  responsible  for  this  anomalous  behaviorc  it  appears  that 
peer  hording  to  the  EGB  core  may  be  responsible  for  tne  seemingly  low 
values  :>f  failing  loading. 
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FIGURE  2.  SANQWiCH  CYLINDER 
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FIGURE  6.  AXIAL- COMPRESSION  BUCXL/NG  COEFFICIENT 
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FIGURE  12.  CHARACTERISTIC  BUCKLE  PATTERN: 
!/8"  CCA  CYLINDER  UNDER 


AXIAL  COM  PRES  SION 


FIGURE  /J.  INFLUENCE  OF  FREE  LENGTH  ON  FAILING  LOAD 
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FIGURE  16 . 5 TRESS  —STRAIN  CHARACTERISTICS  OF 
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FIGURE  21.  CORRELATION  OF  THEORY  AND 
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